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Fano resonances and band structure of two dimensional photonic structures 
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We show that the frequency spectrum of two dimensional photonic crystals is strongly influenced 
by Fano resonances which can be excited already in the linear array of dielectric cylinders. To 
support this claim, we calculate the transmission of electromagnetic wave through linear array of 
dielectric cylinders and show that frequencies of observed Fano resonances coincides with position 
of narrow frequency bands found in the spectra of corresponding two-dimensional photonic crystals. 
Split of frequency band or overlap of two bands, observed in the band structure of photonic structures 
are also associated with Fano resonances. 

PACS numbers: 42.70.Qs 


I. INTRODUCTION 

Frequency spectrum of photonic crystals consists of 
large number of continuous bands. This is a direct conse¬ 
quence of spatial periodicity of the permittivity 5 which 
defines the structure. In the limit of infinitesimally small 
variations of the permittivity, all frequency bands could 
be constructed from the dispersion relation of electro¬ 
magnetic wave in homogeneous medium by reduction of 
the momenta q to the first Brillouin zone [T]-[3]. With 
increasing permittivity contrast, frequency gaps open at 
the edges of Brillouin zone [4] . We will call the resulting 
frequency bands periodic {V) bands. 

For higher permittivity contrast, the frequency spec¬ 
trum is more complicated. It contains, besides the V 
bands, also other, usually very narrow (almost disper- 
sionle^ frequency bands. As an example, we show in 
Fig. Qa) the frequency spectrum of two dimensional 
square array of dielectric cylinders embedded in vacuum. 
Only five of 17 displayed frequency bands are V bands. 

The aim of this Paper is to explain physical origin 
these additional frequency bands. We prove that they 
originate from Fano resonances I3IS] which can be ex¬ 
cited in linear chain of dielectric cylinders by incident 
electromagnetic wave. We find numerically the spectrum 
of these resonances and demonstrate how each resonance 
develops into narrow Fano {T) band in corresponding 
two-dimensional photonic crystal. 

The band structure shown in Fig. [^a) was calculated 
for the square array of thick dielectric cylinders. For 
other structures, for instance photonic crystals composed 
from thinner cylinders, mutual coupling of and original 
V band occurs which results in irregularities of the band 
spectrum. Two of them, namely the split of the V band, 
and an overlap of two bands, will be discussed later. 

Fano resonances in photonic structures [7] were studied 
mostly in process of interaction of two dimensional pho¬ 
tonic slabs with incident electromagnetic wave [8] and 
were used for experimental identification of spectra of 
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leaky modes in photonic structures [9]. For the case of 
individual dielectric cylinder, Fano resonances were ob¬ 
served as a result of the coupling of Mie resonant states 
mm with an incident electromagnetic wave in spheri¬ 
cal m and cylindrical dielectric object m- Fano reso¬ 
nances play an important role in design of metamaterials 
[niiis] and influence considerably the transport proper- 
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FIG. 1. (Color online) (a) Frequency spectrum of the square 
array of dielectric cylinders with radius R = 0.4a (a is the spa¬ 
tial periodicity in x and y directions) and permittivity e: = 12 
calculated by the plane wave expansion method [T]. In the 
frequency region a/X < 0.9 we found hve V bands (black full 
lines), hve even (dashed red lines) and seven odd (dot blue 
lines) T bands. Inset shows schematically the structure in 
the xz plane, red arrow indicates the propagation direction of 
the electromagnetic wave, (b) Transmission coefheient of the 
plane wave incident on the linear array of cylinders (The ge¬ 
ometry of the problem is shown in (c)). The incident angle is 
^ = 0 (red line) and d — tt/IOO (blue line). Note the one-to- 
one correspondence between resonant frequencies of excited 
Fano resonances and T bands observed in the frequency spec¬ 
trum of the inhnite 2D structure, (c) Transmission of incident 
electromagnetic wave through linear array of dielectric cylin¬ 
ders. (d) Geometry of the photonic structure composed from 
N rows of dielectric cylinders. [N — 24 throughout this Pa¬ 
per.) Both structures displayed in (c) and (d) are inhnite in 
the X direction. 
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ties of disordered systems m- Here, we are interested 
in Fano resonances which could be excited by plane elec¬ 
tromagnetic wave incident to the linear periodic array 
of dielectric cylinders (Fig. [^c)) as a result of interfer¬ 
ence of leaky guided modes of periodic array of cylin¬ 
ders [4] with incident electromagnetic wave. These res¬ 
onances can be characterized by the resonant frequency 
and lifetime which is inversely proportional to the width 
of the resonant peak. Surprisingly, they are narrower 
than resonances excited in individual dielectric cylinders. 
If the electromagnetic wave propagates across N chains of 
cylinders (Fig. [^d)) then resonances excited in individ¬ 
ual chains couple together and create narrow transmis¬ 
sion band obtained in the frequency spectra of photonic 
crystals. 

II. STRUCTURE AND METHOD 

We concentrate on periodic photonic structures com¬ 
posed from dielectric cylinders parallel to the ^ axis. 
Cylinders possess a frequency independent permittivity 
s = 12 and permeability /i = 1. The embedding medium 
is vacuum. Radius of cylinders is R and the distance 
between two nearest-neighbor cylinders is a. 

The first structure is the linear chain of dielectric cylin¬ 
ders lying in the y = 0 axis (Fig. [^c)). Incident electro¬ 
magnetic wave of wavelength A excites in this structure 
Fano resonances which manifests themselves as sharp ir¬ 
regularities in the frequency dependence of the transmis¬ 
sion coefficient visible in Fig. [^b). Transmission co¬ 
efficient is calculated also for finite photonic slab con¬ 
structed from A = 24 rows of the same cylinders located 
in planes y = (n — l)a (n = 1, 2,..., 24, see Fig. [^d)). 
Observed transmission spectra will be compared with the 
band structure of an infinite square array of cylinders cal¬ 
culated by plane wave expansion method [T]. 

Transmission coefficient is calculated by the transfer 
matrix method [m ED]. For more detailed analysis of 
the spatial distribution of electric field, we use another 
algorithm based on the expansion of electromagnetic field 
into cylinder functions [niiiiHisi- If an incident electro¬ 
magnetic wave is polarized with electric field Ez parallel 
to cylinders (this polarization is considered throughout 
this Paper) then the intensity of electric field scattered 
at a cylinder centered in {xy) = (0,0) can be expressed 
as 

Ef{r,ip)=Jo{r)a^ + 2'^a'^ Jkir) cos{kip) 

k>0 

+ 2iy^a^ Jk (r) sm(kip) 

k>0 QN 

='Ho(r)/3^ + 2j^/3+'Hk(r) cos(k(p) 

k>0 

+ 2iY,p^nk{ r) sm{kip) 
k>0 

for r < R and r > R, respectively. Similar expressions for 
radial and tangential components of magentic field could 


be derived from Maxwell equations [21]. In Eq. Q, 
Jk{r) = J/c(27rnr/A), Rk{r) = Hk{27rr/X)/Hl{27rR/X), 
[24] J/c, Hk are Bessel and Hankel function, A = 27rc/ci; 
determines the wavelength of electromagnetic field in vac¬ 
uum and n = is the index of refraction. 

For another cylinder, centered at x = and y = 
Uya the field Ez is again expressed by Eq. ^ but with 
new set of coefficients a{nx^ny), P{nx,ny) and cylindrical 
coordinates r and (p associated with the center of the 
cylinder. 

In numerical simulations, we calculate coefficients a 
and P from the requirement of continuity of tangential 
components of the intensity of electric and magnetic field 
at the boundary of cylinders. Note that spatial period¬ 
icity of the structure along the x direction considerably 
reduces the number of unknown coefficients since coeffi¬ 
cients a{nx,ny) and P{nx^ny) fulfill the Bloch theorem 

a{nx,ny) = , /3{nx,ny) =/3(0, 

( 2 ) 

where q is the transverse component of the wave vector. 
This enables us to reduce the number of unknown coeffi¬ 
cients f3^{ny) to N X {2Nb + 1), where Nb is the highest 
order of Bessel function used in numerical calculations 
{Nb = 12 in most cases). To find the transmission co¬ 
efficient T, Poynting vector is calculated on the opposite 
side of the structure. Details of the method are given 
elsewhere [25] . 


III. RESULTS 
A. Isolated Fano bands 

Consider first an array of thick dielectric cylinders with 
radius R = 0.4a. Figures [^b) andj^a) shows the trans¬ 
mission coefficient of plane electromagnetic wave propa¬ 
gating through linear chain of cylinders. A series of very 
narrow resonances could be identified. Similar maxima 
and minims has been found numerically in the reflec¬ 
tion coefficient m- Following |8] we interpret these res¬ 
onances as Fano resonances which results from the inter¬ 
ference of incident plane wave with leaky guided modes 
excited in the periodic cylinder row [4|. Indeed, these 
Fano resonances are accompanied by sharp maxima of 
coefficients P defined in Eq. (Eigure |^b)). More de¬ 
tailed analysis of resonances is given in Sect. El 

Comparison of Eiguresj^a) and (b) confirms that Eano 
resonances develop to narrow Eano bands in the spec¬ 
tra of 2D structures. This is also shown in Eigs. ^c,d) 
which present the frequency dependence of the transmis¬ 
sion coefficient of plane wave propagating through the 
slab composed from A = 24 rows of cylinders. 

Different character of V and E bands is clearly visible 
from the spatial symmetry of the electric field Ez shown 
in Eig. The top panel displays the field Ez for frequen¬ 
cies chosen in the center of five lowest V band shown in 
Fig-d c). As expected, the field symmetry changes when 
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FIG. 2. (Color online) (a) The transmission coefficient of 
the linear chain of cylinder with radius R = 0.4a shown in 
Fig. [^c). Only even resonances are excited since the elec¬ 
tromagnetic wave propagates perpendicularly to the cylinder 
chain [261127] . (b) The frequency dependence of parameters 
. (c) Transmission coefficient for iV = 24 rows of cylinders. 
Fano resonances observed in linear chain develop themselves 
to narrow frequency bands, (d) Detail of the the frequency 
dependence of the transmission coefficient shown in (a) and 
(c) in the vicinity of the Fano resonance of coefficient 67 
(a/A = 0.562). 
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FIG. 3. (Golor online) Gomparison of the field distribution 
Ez in dielectric cylinder for frequencies in the center of the 
five lowest V bands (top) and five lowest T bands (bottom). 
Electromagnetic wave propagates through photonic slab of 
the thickness of 24 rows of cylinders. For each frequency, only 
field inside the 10th cylinder along the y direction is shown. 


B. Overlap of two bands 

For thinner cylinders, the V band and the T band can 
overlap. Then, the resulting frequency spectrum depends 
on the g'-dependence of the frequency in two bands and on 
the strength of their mutual coupling. Consider a simple 
model of two bands 2Vi cosg and 2V2 cosg, (centered at 
the same frequency for simplicity) coupled together with 
coupling constant 2^4. Resulting spectrum has a form 

wi,2(g) = (v+r2) cosg±vWTcw^FFc^ (3) 

Two bands cji(g) and cj2(^) are separated by gap if 
4V1V2 < \A\‘^. This happens either when \A\ is large 
or when one of two bands is narrow. 

Band splitting. We found the above mentioned band 
splitting in the frequency spectrum of the square array 
of thin {R = 0.1a) dielectric cylinders (Fig. [^. The 
transmission coefficient through a linear chain of cylin¬ 
ders (Fig. Bb)) decreases to zero for the frequency 
a/A ~ 0.67. This decrease is accompanied by an increase 
of coefficient Pq (Eq. [^. We interpret this decrease of 
the transmission as a result of excitation of broad Fano 
resonance. 

Figure l^c) shows the band structure of the square ar¬ 
ray of thin cylinders. Note that both the second and the 



the frequency increases from one frequency band to the 
next one ma. The bottom panel shows the field 
for the five lowest resonant frequencies identified in Fig. 
[^b). As shown in Fig. |^d), these frequencies corre¬ 
spond to the center of T bands. The symmetry of the 
field is unambiguously determined by the order of excited 
Fano resonance. 


FIG. 4. (Golor online) (a) Linear chain of thin dielectric cylin¬ 
ders (R — 0.1a) and (b) Transmission coefficient T and the co¬ 
efficient /3o as a function of the frequency a/A. Owing to Fano 
resonance the transmission decreases to zero at a/A ~ 0.67. 
(c) The band structure of two-dimensional square array. Fano 
resonance causes the broad gap between the second and the 
third bands, (d) The transmission coefficient of the electro¬ 
magnetic wave propagating through an array of 24 rows of 
cylinders, (e) The intensity of electric field Ez inside six cylin¬ 
ders along the y direction at band edges X 2,3 and F 2 , 3 - 
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FIG. 5. (Color online) (a) Frequency spectrum of a square 
array of dielectric cylinders with radius R = 0.3a. Solid lines 
represent bands with even spatial symmetry, dot lines cor¬ 
respond, similarly to Fig. [^a), to odd Fano bands. Note 
an overlap of the 3rd and 4th bands, shown in detail in the 
right panel. This overlap results from the coupling of two 
bands displayed by green dot-dashed and red dashed line, (b) 
Frequency dependence of the transmission coefficient for the 
linear chain of cylinders. Three Fano resonances could be eas¬ 
ily identified, (c) Top panel shows the transmission coefficient 
for 24 rows of cylinders as a function of frequency a/A. Two 
lower panels present detail of the frequency dependence of the 
transmission coefficient in the vicinity of the Fano resonance. 

m 


third bands have a minimum at the X point. These two 
bands result from the coupling of the V band with broad 
T band displayed by dot dashed (X 2 F 3 ) and dashed 
(F 2 X 3 ) lines, respectively. 

This statement is supported also by the analysis of 
the spatial symmetry of electric field within two bands 
shown in Fig. [^e). Note that the symmetry of electric 
field changes along the line X 2 F 2 . Also, the field in points 
X 2 and F 3 have the same symmetry; of course, the same 
holds for pair X 3 and F 2 . While the field at X 2 and F 3 
possesses the symmetry of the V band, field close to inner 
band edges X 3 and F 2 has a symmetry of excited Fano 
resonance. 


Band overlap. The overlap of V and T bands is ob¬ 
served in the band structure of the infinite square array 
of cylinders with radius R = 0.3a displayed in Fig. Uta) 
As shown in the right panel, the overlap of the 3rd and 
4th bands can be interpreted as a result of coupling of 
two bands shown by dot dashed and dashed lines, re¬ 
spectively. This overlap can be identified also from the 
complicated frequency dependence of the transmission 
coefficient [28] shown in Fig. [^b) for the slab of = 24 
rows of cylinders. Similar analysis could be done for the 
overlap of the 5th and 6 th bands in Fig. Ha). [ 22 ] 


IV. FANO RESONANCE 


In previous Section, we have shown that Fano reso¬ 
nances excited in linear array of dielectric cylinders cre¬ 
ate the T bands in spectra of photonic crystals. Now we 
will discuss physical origin of Fano resonances. 

Fano resonances have been observed recently in the 
most simple dielectric structure - the single dielectric 
cylinder m- If the frequency of incident electromagnetic 
wave coincides with the eigenfrequency of any cylinder 
leaky eigenmode, the last can be excited. Then, the elec¬ 
tromagnetic field in the neighbor of the cylinder is given 
by a superposition of two fields with the same frequency: 
the incident plane wave and field radiated by excited res¬ 
onance m- The excitation of resonance manifests itself 
as a maximum of coefficient (3 shown in Fig. 6 (a). The 
width of the resonance is proportional to inverse of its 
lifetime. Similarly, excitation of resonant guided mode 
in linear array of cylinders can be identified from sharp 
maxima in frequency dependence of corresponding coeffi¬ 
cient (3 (Figs. 2(b), 6 (b)). The interference of two modes 
is the responsible for narrow maxima and minima in the 
transmission coefficient of incident electromagnetic wave 
displayed for instance in Figs. 1(b) and 2(a). 

We start with the comparison of the frequency depen¬ 
dence of coefficients /3 for isolated dielectric cylinder (Fig. 
6 ffi) ) and for an infinite periodic array of cylinders (Figs. 
2[h) andj^b)). For single cylinder and Ez polarization, 
Pk is a solution of system of linear equations m 


f Jk —3-Lk{R) \ ( Oik 
V J'k -cjypk 



( 4 ) 


where Jk = Jk{‘2.TiR/X)^ ( = y/ja/e is an impedance, and 
the r.h.s is given by the expansion of incident plane wave 
into Bessel functions m 


gi2,ri?/Asin0 = ^ (27ri?/A)e*'='^ 

k 


( 5 ) 


As shown in Fig. [^a), Fano resonances of cylinder lye 
very close to those of an array of cylinders (Figs. 2(b) and 
Fig. [^b)). The first two resonances {k = 0 and k = 1 
are relatively broad, especially for an infinite number of 
cylinders. 

On the other hand, Fano resonances with k > 2 are 
significantly narrower when excited in an infinite linear 
chain of cylinders than in individual cylinder. To explore 
how the shape of the resonance depends on the num¬ 
ber of cylinders, we analyze the scattering of incident 
electromagnetic wave on finite cluster consisting from M 
cylinders along the x direction (Fig. 1 (c)). Coefficients 
were calculated as a solution of system of linear 
equations 


Ap = a ( 6 ) 

where A is a matrix of the size M x {2Nb +1) and vector 
a represents incident electromagnetic wave. Figure [^c) 
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FIG. 6. (Color online) (a) Coefficients [3 for (a) single di¬ 
electric cylinder with radius R = 0.4a. (b) linear row of 

cylinders, (c) Coefficients (red) and an inverse determi¬ 
nant (black) for the clusters of M dielectric cylinders 

(M = 1, 2,4, 24,48, oo) along the x axis. Thick blue line 
shows the position of the T band 0.4067 < ajX < 0.4254 
displayed in Fig. 1(a) by red dashed line. Multiple points 
for a given frequency corresponds with different intensity of 
excited resonance at individual cylinders as shown in (d) for 
the cluster with M = 48 cylinders. 

presents coefficients P 2 finite cluster of M cylinders 
excited by perpendicularly incident plane wave. The res¬ 
onance indeed becomes narrower when number of cylin¬ 
ders increases. Since the spatial distribution of electro¬ 
magnetic field along the finite cluster is not homogeneous, 
P 2 acquires various values for individual cylinders inside 
the cluster. As an example, we plot in Fig. [^d) spatial 
distribution of as well as its real and imaginary part 
calculated in cluster of 48 cylinders. 

Figure l^c) shows also the frequency dependence of an 
inverse of the determinant of the matrix A (Eq. which 
determines the eigenfrequencies and lifetimes of leaky 
guided modes [31] exited in the cluster. Comparison of 
the frequency dependence of an inverse determinant and 
P 2 confirms Fano character of observed excitation and 
enables to estimate the sign of the Fano parameter g < 0 
[i|32|. 


V. ABSORPTION 

Finally, we note an another difference between V and 
3F bands: we expect that the 3F bands are much more 


(a) 




FIG. 7. (Color online) (a) Absorption of the electromagnetic 
wave in the photonic structure discussed in Fig. [^but with 
small imaginary part of the cylinder permittivity. Absorption 
is large at band edges, where the group velocity is small, and 
in the region of the Fano bands, (b) Detailed frequency de¬ 
pendence of the absorption in the resonant frequency region 
for Imag e = 0.001 and 0.01. 


sensitive to the absorption loses than the V bands. One 
reason is that typical 3F band is narrow, therefore the 
group velocity of transmitted wave is small. However, 
more important is that T bands are associated with the 
resonance which led to higher intensity of propagating 
electric field. Figure presents the absorption of elec¬ 
tromagnetic field in the array of = 0.3a cylinders with 
small imaginary part of the permittivity. Fano bands 
could be identified from the position of large maxima of 
the absorption. Note that very small imaginary part of 
the permittivity (Imag e/Real e = S.3 x 10“^) causes 
the absorption of 20% of energy when wave propagates 
through an array of 24 rows of cylinders. 


VI. CONCLUSION 

In conclusion, we showed that the band structure of 
square arrays of cylinders can be completely described 
in terms of two kinds of frequency bands. The V bands 
originates from the reduction of the dispersion relation 
to the first reduced zone. The T bands have origin in 
Fano resonances observed in the linear chain of cylin¬ 
ders. The field distribution in R" bands is determined by 
the symmetry of Fano resonance. Two bands, V and R" 
might overlap, which complicates resulting band struc¬ 
ture. Since resonance is typically accompanied by strong 
electric field, we expect that absorption is stronger in R" 
bands than in the V bands. 
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